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Abstract 

Recent inelastic neutron scattering experiments on (VO)2P207 (VOPO) have 
seen previously unanticipated sharp peaks in the dynamic structure factor in 
addition to the pair of triplet modes observed earlier. We argue that the 
additional features are in essence 'shadows' of the previously studied features 
arising due to umklapp scattering, and suggest experimental tests of this pro- 
posal. The basic point is illustrated by some elementary calculations within 
a strong-coupling expansion for the alternating chain with the right geome- 
try (on which the matrix element for umklapp scattering depends sensitively) 
taken fully into account. 

I. INTRODUCTION 

In recent years, there has been a great deal of experimental and theoretical activity 
studying the magnetic properties of insulating compounds that consist of arrays of well- 
isolated one-dimensional magnetic sub-structures. An example is the compound (VO)2P207 
(VOPO), which was initially thought to be an excellent candidate for a spin- ladder com- 
pound based on early experimental and theoretical work Inelastic neutron scattering 
experiments on single crystal arrays established that these early ideas were incorrect and 
the structure consists, to a good approximation, of an array of alternating antiferromagnetic 
chains that are weakly coupled to each other in one direction perpendicular to the chain 
axis. The spectrum of magnetic excitations was mapped out by these inelastic neutron 
scattering experiments. The lowest lying excitation seen is a triplet mode separated by a 
gap from the singlet ground state of the system. This is identified with the basic single par- 
ticle excitation expected theoretically in an alternating chain (see for instance Ref ^ and 
references therein). The experiment also saw an additional triplet mode above this band in 
a large part of the Brillouin zone. The origin of the second mode was unclear, and there 
was some speculation that it could be ascribed to a triplet bound state of the elementary 
excitations that is also expected to exist in these systems 0. Such a state may be stabilized 
by frustrating interactions |^ It is more likely that the two inequivalent magnetic chains 
in VOPO have differing gap energies the second mode would then be the basic triplet 
mode of the second set of chains. 

New inelastic neutron scattering experiments on a single crystal of (VO)2P207 
(VOPO) [0 have been able to map out the dispersion of the basic excitations of the system 
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in greater detail. These experiments however also see additional sharp low-energy modes 
with dispersions different from the modes seen previously. These extra modes are, at first 
sight, extremely surprising and it is tempting to take them to be a signal of some new, and 
hitherto unanticipated features in the spectrum of the system (possibly arising from frus- 
trated couplings between alternating chains). However, we argue that the real explanation 
for the new modes is quite simple: both modes arise from a purely geometric effect having 
to do with the actual positions of the vanadium ions in the unit cell. The two new modes 
may be thought of as shadows of the basic triplet modes arising from umklapp scattering. 

We begin by detailing the geometry involved and use a very simple general argument 
to calculate the matrix element for the umklapp scattering process that is responsible for 
producing a shadow of the basic triplet mode. We then suggest a straightforward check of 
this explanation based on a comparison of the experimentally observed intensities of the 
basic mode and its shadow at various values of the momentum transfer. It is important 
to emphasize at this stage that this check is quite independent of any theoretical estimates 
of the intensity of the basic mode as a function of momentum transfer and relies only on 
relations between experimentally observed intensity ratios; the calculation of the intensity 
of the basic single particle triplet mode as a function of k {2'Kk/h = ki, is the momentum 
transfer along the chain direction, where h is the unit cell dimension along the chain axis, 
which is conventionally labeled the b axis) is a separate problem that has been addressed 
earlier for the simple alternating chain [0] (these results may also be used in conjunction with 
our analysis to give approximate intensities of the shadow of the single particle band, but we 
do not perform that exercise here). A simple consequence of this scenario is the prediction 
that the shadow band will disappear for h = (here 2nh/a = ka is the momentum transfer 
along the crystallographic a axis perpendicular to the alternating chain axis). 

One also expects that shadows of any bound-state mode will also be formed by an analo- 
gous mechanism involving umklapp scattering. We illustrate this by an explicit calculation, 
to leading order in a strong-coupling expansion, of the bound state contribution to the dy- 
namic structure factor for an alternating chain with the right geometry taken into account. 
The calculated intensity ratios do provide an explicit example of the general argument for 
the strength of the umklapp contribution. 

We also briefly explore the possibility that the magnetic interactions felt by even and odd 
dimers (pairs of spins connected by the stronger of the two antiferromagnetic interactions 
in an alternating chain model) are slightly different. We expect that this will change the 
strength of the shadow bands in a significant way. To get a feel for what to expect, we do a 
simple calculation, again within a strong coupling expansion, of the contribution of the basic 
triplet mode to the dynamic structure factor for an alternating chain with the right geometry 
and the small difference in magnetic interactions felt by even and odd dimers. We see that 
this change in the magnetic environments of even and odd dimers leads to a weak intensity 
for the shadow mode even ai h = (in contrast to our result for the simpler alternating chain 
of Fig |l]) as well as a small splitting between the basic mode and its shadow atA;=l/2,3/2 
in the fundamental Brillouin zone ai h = 0. This shadow at h = 0, as well as the splitting 
at k = 1/2, 3/2 are a sensitive test of the difference in the magnetic environments of even 
and odd dimers in the chain. All experiments to date are consistent with the absence 

of a shadow mode at h = 0. However, in the absence of any straightforward symmetry 
reason forcing the even and odd dimers to be equivalent, the possibility that more refined 
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FIG. 1. Geometry of the alternating chain. Note that the staggering of the dimers in the a 
direction is greatly exaggerated in the figure; v « 0.09a, while w « 0.316 with a ^ 7.7 angstroms 
and b ~ 16.6 angstroms. 



experiments with better statistics will see a weak shadow is still open. 

Lastly, it must be emphasized that our entire approach here ignores the frustrating cou- 
plings between chains that has been argued to exist based on the small, but experimentally 
detectable dispersion seen as a function of h. These couplings along the a direction are 
important ingredients of any quantitatively accurate calculation of the expected neutron 
scattering intensity, but are not expected to change significantly any of our conclusions 
regarding the shadow modes. 



II. SHADOW BANDS DUE TO UMKLAPP SCATTERING. 



We begin with a brief review of some of the relevant details of the structure of VOPO: 
Previous work on VOPO has established that the compound may be thought of as an array 
of alternating antiferromagnetic chains (with the chains oriented along the crystallographic 
b axis and the V^"*" ions forming the basic spin 1/2 constituents of the chains) that are 
weakly coupled to each other in the a direction, and essentially decoupled in the c direction. 
As mentioned in the introduction, we will, for the most part, ignore the weak interchain 
coupling as it is not expected to materially change any of our conclusions. Each unit cell 
of VOPO contains eight V^"*" ions, comprising four dimers. Each of the dimers belongs to a 
different chain. There is a pair of chains near z = 0, and a second pair near z = c/2. The 
members of each pair are related to each other via a screw axis transformation and therefore 
there are at most two magnetically distinct chains displaced from each other along the c 
axis. We can thus focus attention on one representative from each pair. These two chains 
have similar (but not identical) structures and the geometrical effects we discuss here are 
very nearly identical for each type of chain. Using the known structure, we can draw out 
the actual positions of the vanadium ions in the a — b plane. To a good approximation. 
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this gives us an alternating chain along the b axis in which successive dimers are slightly 
staggered along the a axis as shown in Fig |l|. There is an extremely tiny displacement of 
the ions relative to each other in the c direction; this is small enough that we feel justified 
in ignoring it in our analysis. Similarly we can ignore tilts of the dimer units away from the 
b axis. 

The magnetic response of a single chain may be modeled (modulo the possible complica- 
tions that form the subject matter of section |rV|) by the simple alternating chain Hamiltonian 

n = ■ S//(z) + gSjiit) ■ Sj{t + 1)] , (1) 

i 

where J is the overall energy scale fixed by the microscopic exchange constants in the system, 
g represents the ratio of the weak and the strong bonds of the alternating chain, and the 
spins are labelled as in Fig |l[ Note that this Hamiltonian is invariant under translations of 
6/2 along the b axis. However, a glance at Fig |1] shows that the staggering of the positions 
of even and odd dimers in the a direction reduces the actual symmetry of the full structure 
to translations by b and not b/2 along the b axis. This of course implies that momentum 
conservation may be violated during a neutron scattering event by integer multiples of 27i/b 
along the chain axis. Note that this is less stringent than the more usual condition (which 
would be in force in the absence of any staggering along the a axis of even and odd dimers) 
that momentum be conserved modulo integer multiples of An/b, and we believe that this 
simple fact is at the root of the observed shadow bands. Thus, we expect that the extra 
modes seen should be displaced by precisely 27r/b from the basic modes of the alternating 
chain. This seems to be the case with the experimental data 0. To clinch the identification, 
we need to be able to make predictions for the intensities of the extra modes relative to the 
basic modes and see how these compare with the experimental numbers for the intensity 
ratios. This is what we turn to next. 

Let us begin our analysis by writing down the usual spectral representation for the 
dynamic structure factor of our system at T = 0: 

SUKuj) = Y.Si^ - Er, + Eo)\{<!>N\S'i-k)\%)f , (2) 

N 

where |$o) is the exact ground state of the system, \^n) is an exact excited state labeled by 
the index A^, Eq and are the energies of the ground and the excited state respectively 
and S^{—k) is defined as 

^^(-k) = (A)i/2^5i(jy'^-^-, (3) 

where the subscript A takes on values I and //, j refers to the dimer index, L is the length 
of the chain and x^-^ is the position of the spin labeled by j and A (see Fig |I]) (here and 
in the rest of our discussion, we will exploit the rotational invariance in spin space to focus 
only on the zz component of the dynamic structure factor). 

It is convenient to formulate our analysis in terms of operators that directly make ref- 
erence to the states of each strongly coupled dimer in the system. This is achieved by 
transforming to the so-called 'dimer boson' representation [p!0|-[T^. Following Ref we 
write the spin operators as: 
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S?{j) = \ [s\j)Uj) + 4(j>(j) - te^p,t'^^{j)t,{j)) , (4) 

S?iU) = \ {-sHj)tM) - tiij)sij) - te^p,tlij)t,ij)) , (5) 

where a, P, and 7 are vector indices taking the values x,y,z, repeated indices are summed 
over, and e is the totally antisymmetric tensor. s"''(j) and t^O) are respectively creation 
operators for singlet and triplet bosons at 'site' j (in the bosonic language, each strongly 
coupled dimer is thought of as a single site; the separation of adjacent sites along the b axis 
is then 6/2). The restriction that physical states of a dimer are either singlets or triplets 
leads to the following constraint on the boson occupation numbers at each site: 



The spin density is given by 



It is also convenient to define 



Note that as the constraint fixes the number of singlet particles uniquely given the triplet 
occupation number, we may as well think only in terms of triplet occupation numbers; we 
will thus refer to any site which is occupied by a singlet as being in the vacuum state. 

Finally, it is useful to note that the alternating chain is readily analyzed in the limit of 
strong alternation (the so called strong coupling limit, with g <^ 1) 0,0. In the language we 
are using here, the lowest lying excitations in this limit are single particle modes (with one 
bare triplet particle excited above the ground state, which may be thought of as vacuum). 
While corrections are certainly introduced to this picture at higher orders in ^f, it is still 
legitimate to think of the basic triplet mode seen in the real system as arising from the 
contribution of the fully renormalized single particle excitation in the system (for a careful 
analysis of this point for the closely related problem of a spin-ladder, see Ref ||13|| ) 



With these preliminaries out of the way, let us now go through the extremely elementary 
general argument for the strength of the umklapp matrix element responsible for the shadow 
bands. We formulate this here only for contributions to the spectral sum coming from the 
fully renormalized single particle states of the system. The basic argument is nevertheless 
expected to remain valid when applied to the contributions coming from two-particle bound 
states; we will see this expectation realized in an explicit calculation later. 

Let us begin by writing the contribution of the fully renormalized single particle states 

as 

Slnk,u;)=J2Si^-^,)MS%-k)\%)f , (6) 

where \q) is the exact, fully renormalized single particle state of momentum q in the chain 
direction (the state of course has z component of its spin equal to 0; we will not be very 



5 



careful in this section about including this information in our notation) and Sq is the energy 
of this state (with the ground state energy set to zero). We may write this state quite 
generally as 

l^) = (^)'/'Ee^'^^'l^iO-)), (7) 

i 

where is the h component of the position vector of the j*^ 'site' (center of the j*^ dimer) 
in the chain and the notation |\l/i(j)) is intended to denote a state that differs from |$o) 
only locally in the vicinity of site j. Furthermore, we can write S'^(— k) as 

^'^"""^ = ( A)i/2(e«/c„./2 ^ e'^^^'0,,{3) + e-^'^"^/^ ^ e'^^'^'Oum , (8) 

j odd j even 

where the operator Ok^{j) is defined as 

OkSj) = ^(cos(^)a.(j) - .sin(^)0,(j)) , (9) 

with w equal to the distance along the b axis between the two spins that form each dimer. 

This now allows us to write the following expression for the matrix element appearing in 
the spectral sum (|^): 

\{q\S%-k)\^o)\' = M|M£|e-'=W2 + e.(fc.-.)V2g.fcW2|2^^^^^ ^ (10) 
where 5 is defined as 



CO 



— X! ^q,kb+2Trn/b , (H) 



n=—oo 



and Ai is given as 

M{q,h) =j:^"'''H'^iij)\OkM\^o) . (12) 
j 

Note that we can make the dependence of Ai explicit by rewriting this as 

M{q, h) = cos{^)MM) + M^)MM , (13) 



where 



= ^Ee"^'^^^^i(j)l0.(o)l$o) . (14) 

Furthermore, previous work has demonstrated that M-^j is identically zero for the single 
particle states of the simple alternating chain Hamiltonian (|1|). 
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Using all of this we can write the contribution of the fundamental single-particle mode 
to the dynamic structure factor as 

Sll'^-iKu) = \M,ih)\'sm\^)cos'i^)6iuj-e,,) h G (0,^) , 

= \M^(k,-A7r/b)\'sm'i^)cos\^)6{u-e,^^,^/,) h e , y ) , (15) 
while the shadow contribution reads 
^ip,shadow(k,u;) = \M^(k, + 27,/b)\'sm\^)sm'{^)Sico - e,,^,^/,) h e (0,^) , 
= \M^ih-2n/b)\'sm'i^)sm\^)6{u;-e,^^,^/,) h G (y , y ) , 

... /T\i2 • 2 / \ . 2 / \ e/ \7 ^GttStT, 

= \M4,{kb- 6n/b)\ sm (— )sm {—)d{uj ~ Ek^^e^/b) h G 

(16) 

Thus, we see quite generally that the intensity of the shadow band should vanish as 
ka (modulo the complications discussed in section Moreover, it is apparent from 
these expressions that the intensity of the shadow is completely determined by the intensity 
of the basic mode as a function of kh. While there are, in principle, a number of ways in which 
this may be checked against the experimental data of Ref §], it is probably best to simply use 
the experimentally observed intensity of the fundamental mode at fixed ka (for values of at 
which the 'dimer coherence factor' sin^(/cfcw/2) is large) to predict the intensity of the shadow 
at the same ka (this avoids complications due to the weak two-dimensional couplings between 
chains that will introduce additional dependence of the intensity on ka). This prediction 
can then be directly tested against the observed intensity of the shadow after correcting 
for effects of the magnetic form factor of the V^^ ion. Note that this procedure makes 
no assumptions about the form of M.if,{k}j) for the alternating chain Hamiltonian (|I|) and 
serves to separate the purely geometric effect leading to the shadow from our approximate 
knowledge of this function. 

Let us conclude this section by noting that entirely analogous arguments can be used to 
relate the expected intensity of the shadow of a bound-state mode to the intensity of the 
bound-state itself (of course, the analog of A^o- is no longer identically zero, but this merely 
complicates the algebra a little). Instead of going through the corresponding argument for 
the bound state modes in detail, we choose to highlight the minor differences involved by 
doing an approximate calculation of the intensity and position of both the bound-state mode 
and its shadow to leading order in a strong-coupling expansion. This is what we turn to in 
the next section. 



III. BOUND-STATE CONTRIBUTIONS WITHIN THE STRONG-COUPLING 

EXPANSION 

The first order of business is to work out the position in the Brillouin zone and the 
energy of the S = 1 bound state formed from the physical (fully renormalized) triplet 
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particles that are the elementary excitations of the alternating chain (|l]). To leading order 
in the strong coupling expansion, this is particularly simple as the physical single-particle 
excitation coincides with the bare triplet particle created by the triplet boson operator as 
far as the energy levels are concerned. Following the approach used in Ref fl^ , it is easy to 
see [0] at leading order that the triplet bound state exists over two separate intervals for 
the center of mass momentum gem (note that the center of mass momentum takes on values 
in the range (0,87r/6)): the first being (47r/36, 27r/6) and the second being (67r/6, 207r/36). 
The energy of the bound state (with the ground state energy set to zero) is given to leading 
order as eB{qcm)/J = 2 — (7(4 cos^(gcm&/4) + l)/4 (these results were first obtained for a 
slightly more general Hamiltonian by Uhrig and Schulz |^). 

The next thing we need is the bound state wavefunction and the ground state wavefunc- 
tion correct to first order in g; note that it does not suffice to know these eigenstates to 
leading (zeroth) order in g as it turns out that the bound state contributes to the dynamic 
structure factor only at first order or higher in g. We will write these eigenstates down 
in the basis of (bare) triplet boson occupation numbers and polarizations. An extremely 



elementary calculation |T^ gives us the following ground state, correct to first order in g: 

\^o) = |o) + f E (l(i)[o], (j + i)[o]) - |(j)[-i], (j + - l(j)[+i], (j + , (17) 

where |0) represents the vacuum state for the triplet bosons, the number in the square 
brackets gives the z component of the spin of the triplet boson and the number in the 
parenthesis gives the site occupied by the boson (two such pairs separated by a comma 
naturally denote a two-particle state in the bosonic Fock space). 

The zeroth order normalized bound state labeled by the center-of-mass momentum qcm 
(and z component of spin equal to 0) can be easily calculated [|l4l to be 



kcm[0]) = E (|(jl)[-l], (J2)[+1]) - |(jl)[+l], (J2)[-1])) , (18) 

i2>ji 

where the bound-state wavefunction fg^^ is given as 



fUjuJ.) = (A)i/2(£^_l)i/2e-«l-?.-<le--H.-.V/^ , (19) 



with e-'^^/2 = 2cos(gcm&/4)- 



Now, following the approach used in Ref [|T^], it is quite easy to see that the 0{g) cor 



rections to this will involve states living in the zero, one, three and four boson sectors of the 
Fock-space for the bare triplet bosons, in addition to a 0{g) correction to the component in 



the two-boson sector. It is an elementary exercise ||Tj] to work out all these corrections except 
the one in the two-boson sector (as this involves first working out the effective dynamics of 
the physical particles to one higher order in g). While the correction in the two-boson sector 
can also be calculated without too much difficulty using the methods referred to earlier, we 
do not bother to do this explicitly here as it is quite clear that this correction term plays no 
role in the leading order calculation of the bound state contribution to the dynamic structure 
factor. In fact, for our purposes here, it clearly suffices to work out the correction in the 
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one boson sector of the Fock space as this is the only correction that affects our calculation. 
This may be written down readily [14| as: 



S^%cmm = -lT.Ur.ij, J + 1) (I(J)[0]) + |(J + 1)[0])) , (20) 

j 

where 5''^-*|gcm[0]) is the first order correction term in the one-boson sector. 

With all this in place, we can begin our analysis of the bound state contribution to the 
dynamic structure factor by writing down the following spectral sum: 

S^:{k,co)= J2 Sicu - eBiqcm))\{qcmmS^{-k)\^o)\' • (21) 

<Jcm=0 

Notice the different range of summation for qcm in comparison with Eqn (^ (it is of course 
understood that the sum is performed only over those sub-intervals in qcm that actually 
support the existence of a S" = 1 bound state). We can now use and (|^) and calculate the 
matrix element appearing in the spectral sum to first order in g using the 0{g) wavefunctions 
calculated above. 

While it is certainly possible to use the notation of section y and only quote the perturba- 
tive results for the analogs of Ai^ and A^^, we prefer to put everything together and directly 
present results for the leading contribution to the dynamic structure factor. Naturally, these 
results are not expected to be quantitatively accurate. Rather, they provide us with a non- 
trivial example of the general argument of section |T| at work. The bound state leads to 
the following basic contribution to the dynamic structure factor for G {An /3b, Sn/Sb) and 
again for kb E (167r/36, 207r/36): 

^bs,bas.^k, ^) = f^(l - 4cos^(^)) sin^(|(^ - ^)) cos^(^)5(^ - esih)) . (22) 

The shadow of the bound state gives for kt G (0, 27r/36), kb G (lOvr/Sfe, lAn/Sb) and kb G 
(227r/36, Sn/b): 

^bs,shadow(k,..) = 1^(1 - 4sin2(^)) cos^(|(«; - ^)) sin\^)5ico - e^h)) , (23) 

where e\{kb)/ J = 2 — g{Aii\r?{kbb/ 4^) + l)/4 gives us the position of the shadow band. The 
intensity in the bound state mode and in its shadow is depicted are in Fig |^. 

We thus see that the bound state mode also acquires a 'shadow' as anticipated earlier 
on the basis of the general argument. Of course, as mentioned previously, this strong- 
coupling calculation has very little quantitative significance. A quantitative analysis of any 
experimental data on the bound state mode and it's shadow would instead follow the analog 
of the procedure outlined the end of section H with the analog of M.^ included in the analysis 
and the necessary changes made to allow for the fact that the spectral sum is to be carried 
out over a different range from the single particle case (note that this type of analysis can, 
in principle, distinguish between bound state and single particle triplet modes based on the 
different intensity ratios between the modes and their shadows in the two cases). 
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FIG. 2. Intensity of the bound state mode lb, and the intensity of the shadow Ig to leading 
order in the strong couphng expansion. Note that the value of ka is different for the two; in each 
case it is chosen to maximize the intensity. The intensities are both normalized by the average 
value of the intensity in the single-particle mode at this order in the strong coupling expansion. 
The value of g used is 0.7, which is approximately right for VOPO [Q] 



IV. A MORE COMPLICATED MAGNETIC HAMILTONIAN? 



In this final section we briefly consider a possible complication that will affect our results 
at a qualitative level: namely the possibility that the magnetic interactions felt by the 
even and the odd dimers are slightly different. This will clearly change the intensities and 
the dispersions of the various modes observed as the magnetic Hamiltonian will now be 
invariant under translations of b and not 6/2 along the b axis. Thus, the staggering of 
the dimer positions along the a axis will no longer be the only thing breaking the larger 
symmetry of translations by 6/2. Clearly, in such a situation, we expect that the 'shadow' 
band intensity will be non-zero even at fca = (indeed we expect that it will depend quite 
sensitively on the difference in the magnetic interactions of the even and odd dimers). To 
get a feel for what to expect, let us work out the intensity of the basic single particle mode 
and its shadow to leading order within a strong coupling expansion. 

The Hamiltonian we have in mind can be parameterized as: 

n = jJ2[i^ + 9^)^iij) ■ S//(j) + ^7(1 + /i)S//(j) ■ SjU + 1)] 

j odd 

+ J ^ [(1 - gX)Sj{j) ■ SuU) + 9{l - f^)Sjj{j) ■ Sj{j + 1)] , (24) 

j even 

where we have in mind that A and /i are both small parameters that model the small 
differences in the magnetic properties of the even and the odd dimers. 
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The calculation of the 0{g) single-particle contribution to the dynamic structure factor 
is quite elementary and involves nothing new. We will therefore be correspondingly brief. 

We begin our analysis by noting that it is now more natural to count states somewhat 
differently. We will restrict the momentum carried by the single particle state to lie in the 
range (0, 2iT/b) and allow for two distinct bands of single particle states labeled by + and — 
subscripts. 

The energies of these two bands are easily worked out to leading order to be 

^±{q) = 1 T |(4A2 + cos\qb/2) + fx^ sm\qh/2)f'^ . (25) 

Moreover, it is quite elementary to see that the corresponding eigenstates to leading order 
in the strong coupling expansion are (again we choose to write down the state with 3^ = 
as this is what we need to calculate the zz component of the dynamic structure factor): 

= {^y^'V2V{q){Y: e"^^'m)+y^{q) E e^'^'m)) , (26) 

j odd j even 



where V{q) = 1/y 1 + \y±{q)\'^ and y±{q) is given as 

_ 2A ± (4A' + cos'(g&/2) + /i^ sm\qb/2)Y/^ 
y^^'i'- cos{qb/2) + tfism{qb/2) ^ ^ 

One other thing we need is the ground state for this model, correct to 0{g). This can 



also be worked out quite easily to be [|14 



\%) = io) + E mm, ij + i)[o]) - i(j)[-i], (j + - (j + 

j odd 

+ E mm, u + i)[o]) - i(j)[-i], (j + - (j + 

j even 

(28) 

We can now use all of this to work out the one particle piece of the dynamic structure 
factor. The resulting expressions are quite messy for general k and not particularly illumi- 
nating. We will write them down here only for the special case of /ca = 0, as this is where 
we expect a real qualitative difference due to the complications we have introduced into the 
problem: 



sin {kbw/2) 
4 



2G^cos{hb/2) g , 2Q^Go 



;i + ° o + ^(cos(m/2) + 



X 6(10 -e^ih)) , (29) 

where we have defined as 

Q^(kb) = cos^{hb/2) + /i^ sm^{hb/2) , (30) 

and G as 
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FIG. 3. Position of the two single particle bands when the interactions felt by the even and 
odd dimers are different. As explained in the text, the 'basic mode' and it's 'shadow' are actually 
both hybrids made up of the two single particle bands in the problem. We have chosen, for 
illustrative purposes, the values A = 0.1 and fi = 0.1. The parameter g is set equal to 0.7, which 
is approximately correct for VOPO 

G±{kk) = 2\±^A\^ + Q\h) . (31) 

The details of the above expressions are not particularly important. We only wish to use 
the above to arrive at some general qualitative conclusions about the nature of the expected 
intensity at various points in the Brillouin zone. The first of these is of course that we 
have some non-zero intensity at the shadow positions even at /c^ = 0. In this context, it is 
important to note that both the 'basic mode' and the 'shadow' are actually hybrids made up 
of the + band and the — band. For small enough /i and A, the intensity switches between 
the two in such a manner that we have one approximately continuous basic mode and 
another much weaker shadow mode that is also approximately continuous. These results are 
summarized in Fig || and Fig |[ Of course, the avoided level crossing between the two bands 
leads to a small jump in position of both the basic and the shadow mode situated near Ti/b 
and 371 /b. Thus the intensity at the shadow positions and the gap introduced by the avoided 
level crossing are sensitive indicators of the difference between the magnetic environments of 
even and odd dimers. As mentioned earlier, all experiments to date are consistent with 

the absence of extra modes at ka = 0, but in the absence of any straightforward symmetry 
reason, more experiments at /c^ = with better sensitivity and statistics are necessary before 
one can completely rule out the existence of such complications in the magnetic Hamiltonian 
of the system. 
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FIG. 4. Intensity at /cq = of the two single particle bands when the interactions felt by the 
even and odd dimers are different. For clarity, the sin'^{ki,w/2) modulation of the intensity is not 
included in the plots. The values of A, /x and g are set as in Fig ^. 

V. CONCLUSION 



The calculations presented here show that simple geometric effects can lead to the pres- 
ence of the shadow bands observed recently in VOPO. These modes are similar to the 'optic' 
modes that can arise generally in coupled alternating-chain systems with more than one 
dimer per unit cell, for example the 'chain' layers in Sri4Cu2404i [|l^. However, in VOPO 
the shadow modes arise from a single chain. It should be possible to test our proposal by 
comparing the experimentally observed intensity ratios with our predictions. Moreover, as 
mentioned earlier, the intensity ratios can also distinguish between single-particle and bound 
state modes The possible contribution of the triplet bound state to the spin dynamics in 
VOPO remains an open question at this time. 
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